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Abstract 

A compendium of thirty previously published open problems in computational geometry is presented. 
[SIGACT News, 32(3) Issue, 120 Sep. 2001, 63-72.] 

The computational geometry community has made many advances in the relatively short (quarter-century) 
of the field's existence. Along the way researchers have engaged with a number of problems that have resisted 
solution. We gather here a list of open problems in computational geometry (and closely related disciplines) 
which together have occupied a sizable portion of the community's efforts over the last decade or more. We 
make no claim to comprehensiveness, only that were all these problems to be solved, the field would be greatly 
advanced. All the problems have appeared in earlier publications, but we believe all remain open as stated. 
We present them in condensed form, without always defining every technical term, but in each case providing 
at least one reference for further investigation. Our list consists of predominantly theoretical questions for 
which the problem can be succinctly stated and the measure of success is clear. We do not attempt here 
to list the wealth of important problems in applied and experimental computational geometry now being 
addressed by the community as it responds to the application-driven need for practical geometric algorithms; 
we hope that an ongoing project to compile a more comprehensive list will address this omission. We encourage 
correspondence to correct, extend, and update a Web version of this list.Q 



1. 



Can a minimum weight triangulation of a planar point set — one minimizing the total edge length — 
be found in polynomial time? This problem is one of the few from Garey and Johnson [GJ79 whose 
complexity status remains unknown. The best approximation algorithms achieve a (large) constant times 
the optimal length |LK9(;|; good heuristics are known |DMM95|. If Steiner points are allowed, again 
constant-factor approximations are known |Epp94, CL98|, but it is open to decide even if a minimum- 
weight Steiner triangulation exists (the minimum might be approached only in the limit). 

What is the maximum number of combinatorial changes possible in a Euclidean Voronoi diagram of a 
set of n points each moving along a line at unit speed in two dimensions? The best lower bound known 
is quadratic, and the best upper bound is cubic |^ A95, p. 297]. If the speeds are allowed to differ, the 



upper bound remains esse ntially cubic | AGMR98|] . The general belief is that the complexity should be 
close to quadratic; Chew [ Che97 | showed this to be the case if the underlying metric is Li (or Loo)- 



What is the combinatorial complexity of the Voronoi diagram of a set of lines (or line segments) in three 
dimensions? This problem is closely related to the previous problem, because points moving in the plane 
with constant velocity yield straight-line trajectories in space-time. Again, there is a gap between a lower 
bound of ri{'n?) and an upper bo und that is essentially cubic [^ha94 for the Euclidean case (and yet is 
quadratic for polyhedral metrics [ BSTY98 1 ) . A recent advance shows that the "level sets" of the Vorono i 
diagram of lines, given by the union of a set of cylinders, indeed has near-quadratic complexity ||AS00b |. 

J3? 



4. What is the complexity of the union of "fat" objects in 



The Minkowski sum of polyhedra of n 



vertices has complexity 0(n^+'^) |AS99|, as does the union of n congruent cubes |PSS01]. It is widely 
believed the same should hold true for fat objects, those with a bounded ratio of circumradius to inradius, 
as in does in [ |ESO0[ . 
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10. 



11. 



12. 



13. 



14. 



15. 



16. 



Can the Euclidean minimum spanning tree (MST) of n points in Mf^ be computed in time close to the 
lower boimd of Q{n\ogn) | GKFS96[? An MST of a connected graph can be computed in time nearly 
linear in the number of edges |Cha97|, but this is quadratic in the number of vertices n for geometric 
graphs. And indeed the best upper bounds for the Euclidean MST approach quadratic for large d, 
e.g., [CK95|. This problem is intimately related to the bichromatic closest pair problem [AESW91 . 



What is the complexity of computing a minimum-cost Euclidean matching for 2n points in the plane ? 
An algorithm that achieves the minimum and runs in nearly 0{n^ '^) time has long been availab le | Vail 
Recently Arora showed how to achieve a (1 + e)-approximation in n (log?i)0(i/^) time |Aro98| . 



What is the maximum number of fc-sets? (Equivalently, what is the maximum complexity of a /c-level 
in an arrangement of hyperplanes?) For a given set P of n points, S' C P is a k-set if \S\ — k and 
S = P n H for some open halfspace H. Even for points in two dimensions the problem remains open: 
The maximum number of fc-sets as a function of n and k is known to b e 0(nfc ^/'^) by a recent advance 
of Dey [ Dey98| , while the best lower bound is only slightly super linear |TotOO|. 



Is linear programming strongly polynomial? It is known to be weakly polynomial, exponential in the 
bit complexity of the input data |Kha80, Kar84|. Subexponential time is achievable via a randomized 
algorithm |MSW96|. In any fixed dimension, linear programming can be solved in strongly polynomial 
linear time (linear in the input size) [ Dye84 , Meg84 ] . 



Can every convex polyhedron be cut along its edges and unfolded flat to a single, nonoverlapping, simple 
polygon? The answer is known to be NO for nonconvex polyhedra |BDE+01|, but has been long open 
for convex polyhedra |She75||O'R0C|. 



Is there a deterministic, linear-time polygon triangulation algorithm significantly simpler than that of 
Chazelle [Cha91]? Simple randomized algorithms that are close to linear-time are known |Sei91|, and a 
recent randomized linear-time algorithm |AGROO avoids much of the intricacies of Chazelle's algorithm. 
Relatedly, is there a simple linear-time algorithm for computing a shortest path in a simple polygon, 
without first applying a more complicated triangulation algorithm? 



Can the class of 3-SUM hard problems G095 | be solved in subquadratic time? These problems can be 
reduced from the problem of determining whether, given three sets of integers, A, B, and C with total 
size n, there are elements a £ A, b € B, and c € C such that a + b = c. Many fundamental geometric 
problems fall in this class; e.g., computing the area of the union of n triangles. 

Can a planar convex hull be maintained to support both dynamic updates and queries in logarithmic 
time? Re cently t he O(log^n) barrier was broken with a 0(log^~'''^ n) update- a nd O(l ogn) query-time 
structure Cha99|] , and the update time further improved to O(lognloglogn) in BJOO |. 



Is there an 0(?i)-space data structure that supports 0(log7i)-time point- location queries in a three- 
dimensional subdivision of n faces? Currently 0{n log n) space and 0(log^ n) queries are achievable 1 3no97| 



Is it possible to construct a binary space partition (BSP) for n disjoint line segments in the plane of 
size less than 0(nlogn)? The upper bound of O(nlogn) was established by Paterson and Yao |PY90|. 
Recently Toth jTotOll improved the trivial fl{n) lower bound to r2(nlogrt/loglogn). Can the remaining 
gap be closed? 

What is the best output-sensitive convex hull algorithm for n points in M''? The lower bound is fl(n log f+ 
/) for / facets (th e outpu t size). The best upper bound to date is 0(n log / -I- (nfy^^log"'^^^ n) with 
d = l/([(i/2j + 1) |Cha96|, which is optimal for sufficiently small /. 



Can the number of simple polygonalizations of a set of n points in the plane be computed in polynomial 
time? Certain special cases are known (e.g., monotone simple polygonalizations [ ZSSM96[ |), but the 
general problem remains open. The problem is closely related to that of generating a "random" instance 
of a simple polygon on a given set of vertices. Heuristic methods are known and implemented | AII9(| . 
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17. 



18. 



19. 



20. 



21. 



24. 



25. 



26. 



27. 



Given a visibility graph G and a Hamiltonian circuit C, determine in polynomial time whether there is a 
simple polygon whose vertex visibility graph is G, and whose boundary corresponds to G. The problem 
is not even known to be in NP [ p'R9£ |, although it is for "pseudo-polygon" visibility graphs [ pS97 |. 



When a collection of disks are pushed closer together, so that no distance between two center points 
increases, can the area of their union increase? It seems the answer is NO, but this has only been 



settled in the continuo us- motio n case |BS98 . The corresponding question for intersection area decrease 
is similarly unresolved | Cap96 1 . 



What is the complexity of the vertical decomposition of n surfaces in M'', d > 5? The lower bound 
of r2(n^]_was nearly achieved up to d = 3 ASOOa p. 271], but a gap remained for d > 4. A recent 



result Kol01| covers d = 4 and achieves 0{i 



for general d, leaving a gap for d > 5. 



What is the complexity of computing a spanning tree of a planar point set having minimum stabbing 
number? The stabbing number of a tree T is the maximum number of edges of T intersected by a line. 
Any set of n points in the plane has a spanning tree of stabbing number 0{^/n) | Aga92 , ChaSS, Wel93 |, 
and this bound is tight in the worst case. However, nothing is known about the complexity of computing 
a spanning tree (or triangulation) of minimum stabbing number, exactly or approximately. 

Can shortest paths among n obstacles in the plane, with a total of n vertices, be found in optimal 
0{n + hlogh) time using 0{n) space? The only algorithm that is linear in n in time and space is 
quadratic in h | KMM97 |; 0(7ilogn) time, using 0{n\ogn) space, is known HS99 . In three dimensions, 
the Euclidean shortest path problem among general obstacles is NP-hard, but its complexity remains 
open for some special cases, such as when the obstacles are disjoint unit spheres or axis-aligned boxes; 



see [MitOO 



22. Can a minimum-link path among polygonal obstacles be found in subquadratic time? The best algo- 



rithm known requires essentially quadratic time in the worst case [MRW92 . What is the complexity of 
computing minimum-link paths in three dimensions? 

23. How many 7r-floodlights are always sufficient to illuminate any polygon of n vertices, with at most one 



floodlight placed at each vertex? An a-floodlight is a light of aperture a. It was established in [ ECOUX95 



that for any a < n, there is a polygon that cannot be illuminated with an o-floodlight at each vertex. 
When a — TT, n ~ 2 lights (trivially) suffice. So it is of interest (as noted in | UrrOO |) to learn whether a 
fraction of n lights suffice for 7r-floodlights. 

Can an n-vertex polygonal curve be simplified in time nearly linear in n? The goal is to compute a 
simplification that uses the fewest vertices of the original curve while approximating it according to 
some prescribed error criterion. Quadratic-time algorithms have been known for some time |CC9(| and 
a recent algorithm achieves time 0(n^/'^+'^) for a certain error criterion [ AVOO | . In practice, the D ouglas - 
Peucker algorithm is often used as a heuristic; it can be implemented to run in time O(nlogn) [HS94|. 



How efficiently can one compute a polyhedral surface that is an e-approximation of a given triangu- 
lated surface in R^? It is NP-hard to obtain the minimum- facet surface separating two nested convex 
polytopes [DG97|, but polynomial-time approximation algorithms are known (pG95, MS95, AS9^) for 



this case, and for separating two terrain surfaces, achieving factors within 0(1) or 0(log7i) of optimal. 
However, no polynomial-time approximation results are known for general surfaces. 

Given a sufficiently dense sample of points on a surface (technically, an e-sample) , reconstruct a surface 
homeomorphic to the original. This has recently been accomplished for smooth surfaces |ACDLOO|, but 
remains open for surfaces with sharp edges and corners. 

Can the interior of every simply connected polyhedron whose surface is meshed by an even number of 
quadrilaterals be partitioned into a hexahedral mesh compatible with the surface meshing? [BEA+99 



It is known that a topological hexahedral mesh exists [Mit9£, Epp96| | , but despite the availability of 



software that extends quadrilateral surface meshes to hexahedral volume meshes, it is not known if all 
hexahedral cells have planar faces. 
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28. Is the flip graph connected for general-position points in R^^? Given a set of n points in R'^, the flip 
graph has a node for each tetrahedrahzation of the set. Two nodes are connected by an arc if there 
is a 2-to-3 or 3-to-2 "bistellar flip" of tetrahedra between the two simphcial complexes. In the plane, 
the flips correspond to convex quadrilateral diagonal switches; in M'^, a 5- vertex convex polyhedron is 
"flipped" between two of its tetrahedralizations. In the flip graph is connected, providing a basis for 
algorithms to iterate toward the Delaunay triangulation. A decade ago, several [EPW90, Joe91| asked 
whether the same holds in M.^ (when no four points are coplanar), but the question remains unresolved. 

29. Can every convex polytope in M.^ be partitioned into tetrahedra such that the dual graph has a Hamil- 
tonian path? Every co nvex polygon has such a Hamiltonian triangulation^ but not every nonconvex 
polygon does [AHMS96]. The existence of a Hamiltonian path permits faster rendering on a graphics 
screen via pipelining. 



30. We close with Conway's venerable thrackle conjecture, which remains open after more than thirty years. 
A thrackle is a planar drawing of a graph of n vertices by edges which are smooth closed curves between 
vertices, with the condition that any two edges intersect at exactly one point, and have distinct tangents 
there. Conway's conjecture is th at the number edges cannot exceed n. Recently the upper bound was 
lowered from 0{n^/^) to 2n — 3 |LPS95|. Conway offers a reward of $1,000 for a resolution. 



Acknowledgement. We are grateful to Erik Demaine for several helpful comments. 
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